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We consider a microscopic theory of F/S/F trilayers with metallic or insulating ferromagnets. 

^^ , The trilayer with metallic ferromagnets is controlled by the formation of non local pair correlations 

among the two ferromagnets which do not exist with insulating ferromagnets. The difference between 

the insulating and ferromagnetic models can be understood from lowest order diagrams. Metallic 

^ , ferromagnets are controlled by non local pair correlations and the superconducting gap is larger if the 

O ' ferromagnetic electrodes have a parallel spin orientation. Insulating ferromagnets are controlled by 

\^ I pair breaking and the superconducting gap is smaller if the ferromagnetic electrodes have a parallel 

spin orientation. The same behavior is found in the presence of disorder in the microscopic phase 

variables and also in the presence of a partial spin polarization of the ferromagnets. The different 

behaviors of the metallic and insulating trilayers may be probed in experiments. 
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I. INTRODUCTION 



Spin polarized quantum transport has focussed an important interest recently. One of the challenges in this field of 
research is to manipulate correlated pairs of electrons in solid state devices. Several possibihties have been proposed 
recently p-ra] , some of which involve superconductivity and magnetism. 

There is a rich physics occurring at a single F/S interface. For instance it is well established that the superconducting 
order parameter induced in a ferromagnetic metal can oscillate in space [[7|-pi][. In S/F/S Josephson junctions these 
i-rt ', oscillations can induce a change of sign in the Josephson coupling ||ri|. This gives rise to the 7r-state that has 
P5 ' been probed recently in two experiments | |l^,p^ . Another effect taking place at F/S interfaces is the suppression 
O ' of Andreev reflection by spin polarization ]14|]. This has been probed in recent transport experiments, either with 
highly transparent point contacts ||l^ or with intermediate interface transparencies [Q. Other systems such as F/S 
interfaces in diffusive heterostructures have been the subject of several experimental investigations p7|-|2^]. These 
Cn . works have generated many theoretical discussions (see for instance Refs. p^-|30[|). 

J^ ] The specific features associated to transport in multiterminal systems have been discussed recently with various 

-Y-\ ' methods such as Landauer formalism 0, lowest order perturbation for low transparency interfaces B or non per- 
f^*) ' turbative solutions for high transparency interfaces l6J. It was shown theoretically that the conductance associated 
t — « to Andreev reflection is equal to the conductance associated to elastic cotunneling |Q,|j. This could be probed in 
^^ ] future experiments by measuring the conductance as a function of the relative spin orientation of the ferromagnetic 
^T ' electrodes. 

-««.^ ' It is also important to understand equilibrium properties in multiterminal hybrid systems. The proximity effect 

j^ ' at F/S interfaces has been discussed in details recently [^,^. It is well established theoretically that there exists 

Cj ' oscillations of the critical temperature in F/S multilayers as the exchange field and thickness of the F layer are 

I varied pS]. These oscillations of the critical temperature have been probed experimentally in several systems: Nb/Gd 



X3 , multilayers ||,|||, Nb/CuMn multilayers ||§, Nb/Gd/Nb trilayers |37| and Fe/Nb/Fe trilayers |3§. Usadel equations 
^ I have been applied recently to discuss diffusive F/S/F trilayers [Q. F/S/F trilayers have also been discussed recently 
in connection with possible device applications such as a superconducting magnetoresistive memory elements pO[ or a 
superconducting spin switch |^l|. The physics of the F/S/F trilayer with insulating ferromagnets is controlled by pair 
breaking [^2|-Q. Single electron states in the superconductor are coupled to an effective exchange field that cancels 
if the two ferromagnets have an antiparallel spin orientation. As a consequence the superconducting gap is smaller if 
V^ ] the ferromagnets have a parallel spin orientation. 

F/S/F trilayers with metallic ferromagnets have been investigated in a recent work |Q on the basis of effective 
Green's functions. It was found that the physics is not controlled by pair breaking, contrary to the F/S/F trilayer 
with insulating ferromagnets. It was found that with metallic ferromagnets the superconducting gap is larger if 
the ferromagnetic electrodes have a parallel spin orientation 143]. It was proposed that the qualitative physics of 
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multiterminal devices can be characterized by linear superpositions of pair states [ft5[ . We can therefore contrast two 
different situations: 

(i) F/S/F trilayers with insulating ferromagnets are controlled by single electron states. The superconducting order 
parameter is smaller if the ferromagnetic electrodes have a parallel spin orientation. 

(ii) F/S/F trilayers with metallic ferromagnets are controlled by non local pair correlations. The superconducting 
order parameter is smaller if the ferromagnetic electrodes have an antiparallel spin orientation. 

The goal of our article is to discuss F/S/F trilayers with metallic and insulating ferromagnets, as well as a "mixed" 
trilayer with an insulating and a metallic ferromagnet. We use two different approaches, either analytical (with some 
approximations) or based on exact diagonalizations. 

The article is organized as follows. The model is given in section ||, as well as technical preliminaries. Half-metal 



ferromagnets are discussed in section II] . This discussion is extended in section III F to describe an arbitrary exchange 
field. We present exact diagonalizations in section [V. Final remarks are given in section M. 



II. PRELIMINARIES 

A. The model 

We consider throughout the article a superconductor in contact with several ferromagnetic electrodes. The super- 
conductor is three dimensional but a one dimensional geometry will also be used in the numerical simulations. We 
describe the superconductor by a tight binding BCS model in which the electrons can hop between neighboring "sites" 
on a square lattice having a lattice parameter oq. The BCS Hamiltonian takes the form 

HbCS = J2 ^^ {ca,aCl3,a + C^j,^Ca,aj + Yl (^"C+ ^C+ ^ + A^CaaCct) , (1) 

where the summation in the kinetic term is carried out over neighboring pairs of sites. Without loss of generality, 
we assume that the superconductor conduction band is half-filled with therefore kpai^ — 7r/2. We note ep the Fermi 
energy {ep —tiora, half-filled band) and we use also the notation D for the bandwidth. The physics does not depend 
on the details of the band structure. Rather than using a tight-binding model we can also use the free electron 
dispersion relation e(fc) = %^i with ep — ^^^^ the Fermi energy. This dispersion relation is truncated by a high 
energy cut-off e(fcniax) = 21? = 2ep. 

The ferromagnetic electrodes are described by the Stoner model 

Hstoner = JI ^* (^a.^C/?,^ + C+ ^Ca,^ j - /lex ^ yc+,^Ca.^ - C+ ^C^.i j ■ (2) 

(a,/3),cr a 

The case of semi-metal ferromagnets is obtained by considering that the exchange field /icx is larger than the bandwidth 



D. This model with no minority-spin conduction channel is discussed in sections ID and IV. The case of partially 



polarized ferromagnets corresponding to /icx < -D is discussed in section [II F 



In the case of half-metal ferromagnets it will be convenient to use the notation 

t [aakpf 

in which the hopping matrix element is normalized with respect to the Fermi energy. We will use also the notations 
Po for the density of states in the superconductor and p? and p^ for the spin-up and spin-down density of states in 
the ferromagnetic electrodes: 

4 = :^i^ (4) 
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In the case of ferromagnetic metals with two spin channels (see section III F ) we will use the dimensionless parameters 









(7) 
(8) 



B. The method 



We use a Green's function formalism (see for instance 



) to solve the microscopic models given in section [I A 



The first step is to obtain the expression of the advanced and retarded Green's functions G^ ' in terms of the advanced 



A,R 



and retarded Green's functions of the disconnected system g^ ■ . This is done by solving the Dyson equation 



G 



R,A 



9 



R,A 



9 



R,A 



G 



R,A 



(9) 



where the self-energy S contains all couplings of the tunnel Hamiltonian. The Green's functions of the connected 
system incorporate all excursions of the electrons in the ferromagnetic electrodes. The convolution in (0) includes a 
summation over space labels and a convolution of times variables. Since we consider a stationary situation, the latter 
can be transformed into a product by Fourier transform. 

The advanced Green's function takes the following form in the Nambu representation: 



gip{t,t')^-z9{t-t') 



Ca,T(i),C+^(i')|> (K,TW,C/5a(^')}> 

=+,(i),c+,(0}) ({<i(0,c,a(t')}> 



(10) 



where a and /3 are two arbitrary sites in the superconductor. A similar expression holds for the retarded Green's 
function. We adopt the following notation for the Nambu components: 



A,Rr 



[AMr 



The Nambu representation of the density of states pa.pii-o) — -^ 9a ai^) ~ 9a ai^) 



will be noted 



a,/3/ 



\ Pf (^) Pgi'^) 



where p"'^{lu) 
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9a.A^) - 9aA^) and p'}'\uj) 
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Once the advanced and retarded 



Green's functions has been evaluated using (0), we can evaluate the Keldysh component 

G+'- = \i + G" ®€\® g+'- (g) [/ + S ® G-^ 



(11) 
where g^ '~ — 2i7TnF{i^ — Pij)pi,j- The Green's function given by (|ll| ) can be used either to calculate transport 



properties (see for instance |47| ) or to determine the self consistent value of the superconducting order parameter as 
we do in the following (see also ||l[). 

This method can be used to treat non local superconducting pair correlations in the superconductor and in the 
ferromagnetic electrodes (see Ref. M). The pair correlations between two arbitrary sites a and (3 can be characterized 



by the non local Gorkov function 



Gt:pi^) 



The local Gorkov function 



1,2 



G^,>(^) 



can be used to determine 



1,2 



the self consistent value of the superconducting order parameter at any site /? in the superconductor (see |5d,y) via 
the self-consistency condition 
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-U 






(12) 



where U is the microscopic attractive interaction. 



In the following we concentrate on equilibrium properties. Namely, the chemical potentials are identical in all 
electrodes and there is thus no current flow. In this situation the Keldysh Green's function (|l^) simplifies into 

G+q-=nF(a;-/io)(G-^-G«), (13) 

where fiQ is the chemical potential. The calculations based on (n3h will be presented in the main body of the article. 
In Appendix ^ we rederive some of our results by using directly Eq. (pi]). 

The self-consistent value of the superconducting order parameter can be obtained by iterating the process on Fig. H^ 
which starts with a uniform gap profile [A^]. From this gap profile we calculate the propagator ^^^(lj) of the 
superconductor isolated from the ferromagnetic electrodes. From Eq. (^ we obtain the Green's function G^'n. From 
Eqs. (11) or (O) we deduce the Gorkov function [G+'~(w)] , which is used to recalculate the superconducting order 



parameter profile via the self consistency equation (|lj) . 
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FIG. 1. Representation of the successive operations involved in the calculation of the self consistent value of the supercon- 
ducting order parameter. 



C. The different approaches used to determine the self consistent gap profile 

1. Position of the problem 

To discuss F/S/F trilayers we need to find reliable determinations of the self consistent order parameter. It is in 
practice impossible to make an exact analytical calculation of the self consistent order parameter except in the limit 
already considered in Rcf. Ea] where the superconducting gap is uniform in space (the superconductor is smaller 
than the coherence length). The reason why we cannot find exact solutions is the following. Let us start with a 
uniform gap profile A^ = Aq and consider the successive operations on Fig. [I]. Because of the contacts between the 
superconductor and the ferromagnetic electrodes, the Green's functions G^'g are not translational invariant. As a 
result in the next iteration, the superconducting order parameter is not translational invariant. The expression of the 
propagators of an isolated superconductor in the presence of a non uniform superconducting order parameter is not 
known in general. The self consistency relation dT3) is thus a functional relation: 

where the notation [A] means that the right hand side depends on all values of the gap profile. As a consequence, we 
cannot find exact solutions for the gap profile. 

2. Local approach 



We present in section III an approximate analytical treatment in which the functional self consistency relation (hj) 



is replaced by a local relation: 



duj 



^P = -U I ^G;'-''^'{Ap,u^). (15) 



To transform ([14) into (p_5h , we assume that the propagators of the isolated superconductor in the presence of a non 

l,i? , rA,R 



uniform gap has the same energy dependence as the propagators with a uniform gap. The propagators g^ 'g and /^ ' 



depend on an effective gap Aq, 



/3- 



3. Exact diagonalizations 



We present in section [V another possible approach in which we use exact diagonahzations to solve exactly the 
functional form of the self consistency equation (|lj). This numerical method is restricted to small system sizes. We 
will find that the exact diagonalizations are consistent with the "local" approach in section [II in the sense that we 
find Aaf < Ap with both approaches for metallic ferromagnets. 



D. Green's functions in the presence of a uniform superconducting order parameter 



We end-up this preliminary section by giving the form of the Green's function of a superconductor having a uniform 
superconducting order parameter: A^ = Aq for all sites /3. The spectral representation has already been given in 
Ref. 0, as well as the final form of the propagators below the superconducting gap. The final form of the propagators 
above the superconducting gap is found to be 
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where the phase in the prefactor is given by 

lj;(uj) = \Xa- Xp\J{uj - ^isf 

Vf V 

and ip{uj) — kp\xa ~ Xf}\. li lu ^ A, the Green's function reduces to 
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exp[i{ip{uj) +ip{uj))], 



which will be used in section [II D 
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III. HALF-METAL FERROMAGNETS AND FERROMAGNETS WITH BOTH SPIN CHANNELS: 

ANALYTICAL RESULTS 



We discuss in this section ma inly t he solutions of F/S/F trilayers with half-metal ferromagnets having a single spin 
conduction channel. In section III F we extend our discussion to a model having both spin conduction channels. As 
discussed in section II C 1 , the self-consistency equation for the superconducting order parameter (|l4| ) is a functional 
of the gap profile. In this section, we replace the functional equation (14) by the local equation (|l5|). 



A. Superconductor connected to a single-channel half-metal ferromagnet 

We first consider the case on Fig. where a superconductor is connected to a single-channel half-metal ferromagnetic 
electrode. Using the Dyson equation ^ and the expression ( fi3| ) of the equilibrium Gorkov function, we obtain easily 
the local Gorkov function 



G^;p = 2i7:nF{ij) i Pf + ^^ 
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_}_ a,a,A (3.a,Ara,l3.A ^a,a,Rj3,a.Rfa.0,R 



(18) 



with 



25 = 1- |i"'"P.gi;iV'"- 



(19) 
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FIG. 2. Representation of a model in which a single channel spin-up electrode is in contact with a superconductor. 
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FIG. 3. Variation of the superconducting order parameter with a single ferromagnetic channel. We used realistic param- 
eters: the Fermi energy is ep = 10 eV and the value of the attractive electron - electron interaction is such that the bulk 
superconducting order parameter is Abuik ~ 1 meV. 



1. Gap profile with fixed phases 



We first make the additional assumption that the electronic phase in the Green's function 
distance: ip 
be 



61) does not depend on 
-7r/2 for all distances. Phase averaging will be discussed in section HID. The gap profile is found to 
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(20) 



where to is given by Eq. o). Ra^p — \xa — Xjsl is the distance between sites a and /3 in the superconductor. Far 
away from the contact the superconducting order parameter is equal to the bulk value. The minimum value of the 
superconducting order parameter at the contact can be estimated from Eq. ( po| ) by replacing Ra,p by the lattice 
spacing oq: 



2D exp 



U tooq 



[1 + 4] 



The complete gap profile is shown on Fig. ra for several values of the hopping between the superconductor and the 
ferromagnetic channel. 



2. Role of 2kF oscillations 
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FIG. 4. Variation of the superconducting order parameter as a function of the distance to the contact in the case of a 
superconductor connected to a single-channel ferromagnetic electrode. We have incorporated the oscillatory phase factor 
in (|2l|) and we suppose that kpao = ^. The period of the oscillations is thus 2ao, as expected for 2A;_f oscillations. The 
parameters are the same as on Fig. ^ 

Now we discuss the role played by the phase (p appearing in the Green's function g^^p (see Eq.(|l6|)). In the presence 
of this phase factor, the self consistent superconducting order parameter develops 2kp oscillations: 
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(21) 



where Ra.fs = \xa — a?^|- The gap profile is shown on Fig. 0. One may notice that Abuik — ^{Ra,p) deduced from 
Fig. y is related to the wave function of a spin-up electron injected at site a in the superconductor. Namely the 
superconducting gap is maximal when the spin-up wave function is minimal. 



B. Superconductor connected to two single-channel half- metal ferromagnets 
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FIG. 5. The model in which two single channel half-metal electrodes are in contact with a superconductor. 
Now we consider that two single-channel half-metal ferromagnets are connected to a superconductor (see Fi g. |^) 



We assume that the electronic phases are fixed to the value if 
the discussion of phase averaging. 



-7r/2 for all distances and postpone for section 



HID 



1. Antiferromagnetic alignment 

Let us consider the model on Fig. g in which two single-channel half-metal spin-up and spin-down electrodes are in 
contact with a superconductor. The local Green's function takes the form 

The propagators G""'^ and G° '^ are given by 
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and where 
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(23) 

(24) 

(25) 
(26) 

(27) 



The self-consistent superconducting order parameter is obtained by evaluating the high energy behavior of the local 



Gorkov function given by (13). We use a "local" approximation in which the gaps A^jj and Aa'^p appearing in the 
propagators ga.f3 and fa,i3 are replaced by the local gap A^. The high energy behavior of the local Gorkov is found 
to be 
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(28) 



(29) 



where ig'" and ig '" are the tunnel matrix elements normalized to the Fermi energy (see Eq. (^). 



2. Ferromagnetic alignment 



Using the same method for the ferromagnetic alignment, we obtain the high energy behavior of the local Gorkov 
function 
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Comparing Eqs. ([2S|) and (p3|), we see that: 



(30) 



(31) 



(32) 



(33) 



(i) As expected, the "local" contributions of order l/i?^ « and 1/i?^/ a do not depend on the relative spin orientation 



of the ferromagnetic electrodes, 
parameter. 



The local contributions generate a reduction of the superconducting order 



(ii) The lowest order "non local" contribution arises at order l/[Ra,i3Ra',i3Ra.a'], and depends on the relative spin 
orientation of the ferromagnetic electrodes via a factor of two in Eq. ( p3| ) , not present in Eq. ( |29| ) . Because of 
this non local contribution, the superconducting gap is larger in the ferromagnetic alignment, which can receive 
a simple interpretation in terms of the diagrams contributing to this non local term (see section [II E ) . 



3. Gap profiles 
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FIG. 6. Variation of the superconducting order parameter in the presence of two ferromagnetic electrodes. It is assumed that 
the point (3 is aligned with the points a and a' . The coordinate is normalized to the separation Ra^a' between the contacts. 
The different curves correspond to Ra,a' ~ 4 (O), Ra,a' ~ 8 (+) and Ra.a' = 16 (□). We used the same parameters as on 



Fig. M. The contacts have a low transparency: 
orientation. 
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0.1. The ferromagnetic electrodes have an antiparallel spin 
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FIG. 7. The same as Fig. M with high transparency contacts; 



— f 



The gap profiles are shown on Fig. g in the tunnel regime and Fig. [^ in the high transparency regime. The gap 
is reduced close to the contacts with the ferromagnets, which was already obtained for the single channel model in 
sections |IIIA| and |III A 2| (see Figs. | and |). 



4- Difference between the superconducting gaps in the ferromagnetic and antiferromagnetic alignments 



o 

a 

o 

O 

h-1 




-2 2 

Normalized coordinate 



FIG. 8. Variation of the logarithm of S/j defined by (p4J). The ferromagnetic gap is larger than the antiferromagnetic gap. 
The parameters are the same as on Fig. H. The contacts have a high transparency: t^'" = tg '" = 1. 

At each point /? in the superconductor, we define J^g as 
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Sp is positive for metallic ferromagnets, and takes a simple form at large distance: 
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C. Superconductor connected to two single-channel insulating ferromagnets 



Now we show that we recover the correct physics in the case of insulating ferromagnets B3]. The propagator 
relevant to describe a ferromagnetic insulator decays exponentially with distance and is su ch that gf^ = g^^ . The 
local propagators ga,a and ga' a' are real numbers. Using the same method as in section IIIB we find that the Gorkov 
functions are still given by (gsj) and (^0|) but with a different form of A^^ and A^. The expression of the Gorkov 
functions to order 1/ R^ is the following: 
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We deduce from ( p6| ) and (^) the value of (5/5 defined by (^J): 
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The values of the normalized tunnel matrix elements are such that lip'^H^o '" I < 1- As a consequence we recover the 
perturbative result obtained in Rcf. I^^l for insulating ferromagnets ((5/3 < 0). 

D. Phase averaging 

The microscopic Green's function g"'^ depends on the phase variables (fa.pii^) and tpa.i3{u)) (see Eq. (|lq)). In the 
preceding subsections, we have assumed that these phases were fixed to (fa,i3{t^) — — 7r/2 and ipa,i3{^) ~ for all 
distances. In fact the microscopic phases phases given by (fiT[ ) oscillate rapidly on microscopic scales, as opposed to 
the slowly varying prefactor involving 1/ Ra,p in ga,p (see Eq. (|l6|)). Moreover in a multichannel model these phases 
are averaged out when the summation over all channels is carried out (see Ref. Q). It is thus legitimate to consider 
the phases as random variables and to average the Gorkov functions over "phase disorder" . 

1. Single- channel problem 

Let us start with the single-channel problem on Fig. 0. The local Gorkov function were already given in Eq. (fsl) . We 
need to evaluate l^l^gP'°'^^f"-'0'^)) and (^(^gP'°'^^f°''0^^))^ where (( )) denotes the averaging over phase disorder. Assum- 
ing that the phases are symmetric (namely ipa,i3{'~^) = 'Pi3,a{oj) and ipa.,i3{^) = i^f3,a{^)) leads to l^l^gl^^°''^f"'P^^)) — 
and (((^/5'">^/">^'-R)) = 0, from what we deduce that the superconducting gap does not depend on the transparency 
of the contact with the ferromagnet. Since this conclusion is not acceptable physically, we suppose instead that the 
phases are antisymmetric: Lpa,f3 = — 7r/2 + kpRa.p, Vp,a = — 7r/2 — kpRa^p, and 'ipa,i3 = —'4'p,a- Then the expectation 
value of ((c,/3,",^/">/3^^)) is finite, as it should: 
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from what we deduce the self consistent superconducting order parameter 



A/3 = 2D exp { - 



1 27r^h^ 
U ma^ 



ao 



1-' 

2 V i?Q,/3 



1 + io 



(39) 



This form of the gap profile is similar to Eq. (^0|), except for the coefficient 1/2 due to phase averaging. Now we 
discuss phase averaging in different types of two-channel heterostructures. 



2. Superconductor connected to two single-channel half-metal ferromagnets 



With two half-metal ferromagnets, we obtain 

.2 / \^a,y 



A^F - 1 







I \ • 2 / \ ^0 / I I 



Ra.pRa' .pRa.a' \1 
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It 
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I 



K',p Vi + l^o^P 



Sin^ ((<5a',/3) 



(41) 
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ua:"|2 N / UQ .a |2 \ 
R 7? R \ 1 i°i^'^."|2 -; "a' a' 1 9 i^°^ ('^/5.« + fc.,a')sin {ipc',p) + COS {>ff3,o' + >fa',a) sin ((/J^./s)} ■ 

n„_fiRr,'j!n„,y' V 1 + Itg 1/ \ 1 + 1*0 p/ 



After averaging over phase disorder, we find 

((^Mctal)) = 1 ^ oljT" 



\tr\' 



^RIbK^ + K 



1 a2 



|to '" P 



2i?'^^^ I 1 + It)^''"' P 



(42) 
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2 Ra.fjRa' .pRa.a' \ 1 + 1^0 
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1 + I<'"'P. 






(43) 



a,a\2 



\tr\ 



Ra^pRc- ,pRa^c.- Vl+|io IV Vl + I^O " 



The form of the Gorkov functions is therefore similar to Eqs. (|29|) and (p3D, except for the 1/2 prefactors 



5. Superconductor connected to two single- channel insulating ferromagnets 
For a superconductor connected to two insulating ferromagnets, we obtain 



((A^.^» = 1 - ' "° 



2^1/3 



K-y 



i+iirp 



1 a;^ 



,a ,Q 14 



2^«',/3 
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"0 



((AL» = 1 - ' "° 
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2 Ra,pRa' ,l3Ra 

2 / |,a,.|4 



-(i-itrnto''"'i') Y 
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+ 1^: 



a. a I 
I 






^^0 /-|^ _ |^a,Q|2|^a',a'|2\ / I'-Q \ 

Ra,pRa',l3Ra,a' ^ -^ Vl+l^o'"P 

which differs from (Bq) and (|37|) by the 1/2 coefficients. 



to '" P 



l + lil 



I 
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a >" |2 



i + |tS'^"'p 



(44) 



(45) 



^. "Mixed" junction with an insulating and a metallic single-channel ferromagnet 

Now we consider the "mixed" heterostructure on Fig. in which electrode "a" is a single-channel half metal 
ferromagnet and electrode "b" is insulating. Using the same formalism as in the preceding sections, we obtain 



A 



AF 

Mixed 



<, vi + i^r'^ 



K " \' 



Rl'..0 li + l^'"'!' 
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(47) 



\t 



a ,a |2 
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Rl',Ai + K''"'\ 



cos (<^;3,Q') sin {(pa',f3) - 1^0 '" P sin ("l^a',/?) ^in (V5/3,a') 
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^a,f3^a.' ,l3''^a,a.' 



\to'y 



1+ 1] 



I 
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i + K''"'\\ 



[sin {ipi3,a + ipa,a') sin {^a',f3) + sin {ipp.a' + (/Jq'.q) Sill {(pa,0) 

+ |io '°' P (cos ((p;3^Q + 'Pa,a')sVL\ {ipa' ,p) + COS ((/3/3^q' + ((5a', a) sin ((^a./j)) 



Averaging over phase disorder leads to 
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(48) 



(49) 



i?a,/3i?a',/3i?a,a' Vl + ICiV V 1 + l*o''"' I V 

As a consequence this heterostructure behaves hke the full metallic heterostructure {Sp > 0). 

E. Lowest order diagrams 

In this section we point out a simple rule that can be used to determine whether Ajj > A^^ or whether A| < Ag^. 



Let us start with the metallic model. We see from Eq. (p2) and Eqs. (|25|) - (26) that the lowest order non local 
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FIG. 9. Lowest order processes in the case of metallic ferromagnets. (a) corresponds to (pOt), (b) corresponds to (|5l|), and 
(c) corresponds to (p2). 

process are given by 
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a. J.Q ,a 



a ,a la ,q' q .(3 



92.2 



if the ferromagnets have an antiparallcl spin orientation, and by 

g(},a^a.agapa,agaW^a'.a'ga\a'^a',a'ja',l3 



(50) 

(51) 
(52) 



if the ferromagnets have a paraUel spin orientation. The corresponding diagrams are shown on Fig. |9|. Because each 
of these diagrams contains four "g" propagators and one "/" propagator, the sign of the coefficient of the non local 
term in A^^^^g^j and A^^^^^j is positive (in agreement with Eqs. |2^ and |3^). There is one diagram in the case of parallel 
spin orientations, and there are two diagrams in the case of antiparallcl spin orientations, which explains the factor 
of two appearing in the case of a parallel spin orientation - see Eqs. (29) and ( p3| ) 

but now 



Let us now consider insulating ferromagnets. The lowest order diagrams are given by (pOh, (51) and (t 
3°'° and (7° •"" are real numbers. As a consequence the sign of the non local term with insulating ferromagnets is 
opposite to the sign of the non local term with metallic ferromagnets, which is in agreement with Eqs. (pq) and (137 







a. 


Z^ Metal / spin-up 


P V. 




a' C= 


yj Insulator / spin-down 







FIG. 10. Lowest order processes in the case of the mixed junction with antiparallel spin orientations. 

Finally in the mixed case the diagrams given by (pQ), ( |5l| ) and ^2) cancel because they are pure imaginary. 
Therefore we look for the diagrams appearing in the next order. One of these diagrams is represented on Fig. nO. 
There are four "5" propagators involved. The sign of the diagram is positive, which explains why the mixed junction 
behaves like the metallic junction. The diagram on Fig. nffl is proportional to jf^^pji" '" j**, which is in agreement 
with Eqs. (H) and iM). 



F. Ferromagnetic electrodes with both spin channels 



In this section we calculate the superconducting gap of a superconductor connected to one-dimensional ferromagnetic 
electrodes having a partial spin polarization. The motivation is to show that the results obtained in the preceding 
section for half-metal ferromagnets are valid also in the presence of a partial spin polarization. The ferromagnetic 
electrodes are described by the Stoner Hamiltonian (H) in which the exchange field is smaller than the Fermi energy. 
The spin-up and spin-down channels are characterized by the density of states p| and p^ (see Eqs. (pi) and (P)). The 
derivations of the results is given in Appendix p^. 

We use Eqs. (B4) - (B7) and Eqs. (BS) - (Bll) obtained in Appendix M and the local approximation already 
discussed in section pi We find that the high energy behavior of the Gorkov function is given by (^) and ( |30| ) , with 
the following parameters A^|l,'[°j 



and A^F^^j. 
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where x^ and xi are given by Eqs. (|^) and (H). The case of half- metal ferromagnets discussed in section III can be 
recovered by considering the limit x^ = 0, in which case Eqs. ( |5^ ) and (^) reduce to Eqs. (^) and (^). On the 
other hand it is easy to show from Eqs. (^^ and ( |54| ) that the ferromagnetic and antiferromagnetic superconducting 
gaps are equal only if there is no spin polarization (p| — Pi)- As a consequence for two metallic ferromagnets having 
an arbitrary small spin polarization, the ferromagnetic gap is larger than the antiferromagnetic gap. This generalizes 



the behavior obtained in section III in the case of half-metal ferromagnets. 



IV. EXACT DIAGONALIZATIONS FOR HALF-METAL FERROMAGNETS 



We present in this section a simulation based on exact diagonalizations in which we ca n iter ate t he fun ctional form 
of the self consistency e quati on given by Eq. (|l^). The method is presented in sections IV A and IV B . The results 
are discussed in section IV C . 



A. Bogoliubov-de Gennes equations 

1. Boqoliubov-de Gennes Hamiltonian 



L sites 



a a 



a a 



Ferromagnet Superconductor Ferromagnet 

FIG. 11. The geometry treated in the numerical simulation. A superconductor on a one dimensional segment with L sites is 
connected to two ferromagnets. 

We consider the BCS model defined by Eq. (|l|) on a one dimensional lattice with L sites (see Fig. [ll| ): 

L L 



o".i— 1 



i=l 



The two dimensional model cannot be treated numerically because of computational limitations and this is why we 
consider a one dimensional geometry. Nevertheless, the method that we use in ID can be also applied to 2D models. 
It is convenient to use the notation 



i^t 



--1,1' 



=L.T'Cia.--'CL,i 



in which ^t has 2L components. We use (p6[) to obtain the Bogoliubov-de Gennes Hamiltonian 



where the kinetic term is 



H = i^+Ki;^+i^+Ai^^, 



^ki = *(4j+i +4,;-i) 



m4 



K, 



1:2 



K. 



2,1 
kA 



and the pairing term is 



Al,^ - A^f = 



(56) 

(57) 

(58) 
(59) 
(60) 

(61) 
(62) 



Similarly to the Nambu representation, we have used the label "1" for the "electronic" components of -0 and the label 
"2" for the "hole" components. We have doubled the space coordinates: the labels i,j correspond to the electronic 
component and the labels fc, I correspond to the hole component. The symbol 5i_k means that i and k correspond to 
the same site on the lattice but belong to a different Nambu component. 
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2. Spectral representations 
The eigenvectors of the Bogohubov-de Gennes Hamiltonian ( |57| ) take the form 

L L 



i=l fc=l 

L L 

|V'/3)=E^/3^*l^»)+E^'3'^|efe)' (64) 

4=1 fc=l 

where the eigenvalues are such that Aq > and Xp < 0. In this notation there are L kets je;) associated to the first 
component of the Nambu representation, and there are L kets \ek) associated to the second component of the Nambu 
representation. We deduce from Eqs. ( |63| ) and (|6l) the form of the quasiparticle operators 

r+ ^ = ^ Ra,iCi,^ + ^ ^a.feC^i (65) 

i k 

r/3:T = Yl ^/5.»^».T + Yl ^/3.fe4,i (66) 

i k 

which diagonahze the Bogohubov-de Gennes Hamiltonian 

a (3 

The spectral representation of the Green's function ( |lO|) can be expressed in terms of the matrix R: 

A, 1,1 _ V^ Rp,iRfJ.j V^ Ra.iRa,j /„„N 

^''' ~^Lo + ^ll^[^Ji+\Ef,\y^LO + ^l^-[^i-E^] ^ ^ 

A, 1,2 _ V^ Rp,iR(3.k Y^ Ra.iRa,k , „„-. 



B. Evaluation of the Green's functions 

1. Evaluation of a spectral representation: 



The Green's functions are obtained from Eqs. (p7[) and ( pq ) in terms of their poles u;„ and residues i?„: 

5o» = E ^'^^ (69) 

To make the integration over energy, we go to the limit of zero dissipation {rj — > 0) and use the identity l/[a;— a;„— 177] = 
V /[oj — a;„] + iirSiu; — a;„). To show that the principal part can be neglected if w > A, we come back to the particular 
case where the superconducting order parameter is uniform: Ap = Aq for all p. In this case the spectral representation 
takes the form 



k 

We start from Eq. (^Ol) and make the substitution 

1 



a,/3Vw) - 77 >,e "' -— — -H — . (7U) 

^ N^ uj - (lis + Ek) - tri Lu - (us - Ek) - ti] 



nrS {lo ~ {^is + Ek)) (71) 



^ - (/-is + Ek) - iv 

^ . -^zTTS{u;-{fis-Ek)). (72) 

^ - it-is - Ek) - 17] 
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The Green's function given by Eq. (70) becomes 

1 



1^ 



-jfT.^''' 



^>) [{uk)^S [uj - ifis + Ek)) + {vkfS {uj ~ {fis - Ek))] 



(73) 



After using the (5-function (|73| ) and performing the integral over wave vector we recover the form (^^ of the Green's 
function in which the term proportional to cos (p has been discarded. The fact that the term proportional to cos (p is 
not included does not constitute a problem because we know from section [II that the envelope of the 2fcF-oscillations 
is the same in the presence or absence of the cos (p term. Therefore if w > A, Eq. (|69| ) can be replaced by 



30 {u) = iTT 2^ RnS {uj - LOn) ■ 



(74) 



2. Evaluation of the S-functions 



To evaluate the (5-function in (|74|), we replace S{uj — aj„) by 5jj{u}), where (5^(w) is a function having a width r] in 
energy, and normalized to unity: ^ 5.,^{u>)du> = 1. For instance 5,,(w) can be chosen as a Lorentzian or a Gaussian. To 
obtain the value of a Green's function at a single energy w the Lorentzian or Gaussian will be evaluated 2L times. 
To optimize this part of the program, it is useful to use a function (5^(w) that is finite only in the interval [—77, rj\ and 
vanishes outside this energy interval. The simplest choice is given by 



^^M = -^ 



1-1^ 



if 



uj\ <rj. 



(75) 



C. Results 

We consider the geometry represented on Fig. O in which a one dimensional superconductor on an open segment 
with L sites is connected to two ferromagnetic metals. The superconductor is described by the BCS tight-binding 
Hamiltonian dl^). We note <o = ta,a/t = ta',a'/t the tunnel matrix element connecting the ferromagnets and the 
superconductor, normalized to the bandwidth of the superconductor. Low transparency interfaces correspond to 
to <S^ I and high transparency interfaces correspond to io ~ 1- 



1. Density of states 



We have shown on Fig. O the energy dependence of the density of states p"'^ associated to the ordinary propagator 



(see section [IB) in the presence of a uniform gap profile A^ = Aq for all (3, and with L = 128 sites. It is visible on this 
figure that the different parameters of the simulation are compatible with each other. Namely, the level broadening 
77 is sufficiently small to have a well-defined superconducting gap. The level broadening is also sufficiently large for 
quasiparticle states to form a continuous band. Because of these two constraints, we cannot use in this simulation 
realistic parameters as we did for the local approach in section III (see Fig. 0). Using realistic parameters would 
require too large system sizes. 



Finally, the calculation of the superconducting order parameter presented in section III were based on the estimation 



of the Gorkov function at high energy. By contrast low energy degrees of freedom play a relevant role in our simulation. 
One of the questions that will be answered by the exact diagonalizations is to determine whether low energy degrees 
of freedom (probed by the numerical simulation with strong finite size effects) have the same physics as high energy 
degrees of freedom (probed by the local approach in section III). 



2. Gap profile 



32 and L = 
Namely, the 

superconducting order parameter is reduced close to the interface with the ferromagnets and we obtain 2kp oscillations 

in the gap profile. 



The gap profile is shown on Fig. 03 for L = 128 sites. We have obtained similar results for L - 
64 sites. The gap profile obtained with exact diagonalizations is qualitatively similar to section III 
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FIG. 12. Energy dependence of the density of state Pg , for two values of the distance between the sites a and /3. We used 
periodic boundary conditions, with L — 128 sites. The hopping energy is t = 0.5, the superconducting gap Ao = 0.2 is uniform 
and the level broadening is r; = 0.1. 

3. Difference between the superconducting gaps in the ferromagnetic and antiferromagnetic alignments 



We have shown on Fig. |lj the variation of Sp defined by Eq. ( |34| ) with L = 128 sites. Similar results have been 
obtained with L — 32 and L = 64 sites. For each site (3 in the superconductor, we have calculated the superconducting 
order parameters A^^ and A^ with parallel and antiparallel spin orientations in the two ferromagnetic electrodes. 
From what we deduce the value of Sp defined by Eq. (|34|). The calculation consists in iterating the process on Fig. Q until 
a sufficient precision has been obtained. The relative error made in the determination of the order parameters is several 
orders of magnitude smaller than the difference between the ferromagnetic and antiferromagnetic superconducting 
gaps. 

We made the simulations with two values of the normalized hopping between the superconductor and the ferro- 
magnet (see Fig. |lj). We also tried larger values of the interface transparencies but the algorithm did not converge. 
The clarification of this point is left as an open question for future work. From the result presented on Fig. nj we 
deduce that 



(i) With all available sizes and interface transparencies, 6/3 defined by (p4) is positive, meaning that the gap is 
larger in the ferromagnetic ahgnment. This is opposite to the model with insulating ferromagnets M^l and is in 
agreement with the approaches used in sections [11 and [II F . 



(ii) Sp defined by Eq. (p4) tends to zero in the bulk of the superconductor. The cross-over between the surface and 
bulk behavior is controlled by a length scale which is of order 10 on Fig. nj. It is expected that this length scale 
is equal to the superconducting coherence length given by ^o — ^p/ikpAo). 



V. CONCLUSION 



To summarize, we have provided a detailed investigation of F/S/F trilayers with metallic ferromagnets. We found 
that the physics of the metallic trilayer was dominated by pair correlations, not by pair breaking. This behavior was 
obtained with several complementary approaches: 

(i) An approach based on the estimation of the high energy behavior of the Gorkov function for half-metal ferro- 
magnets (see section [II) and for ferromagnets having both spin conduction channels (see section [II F ). In these 
approaches we could use realistic parameters (Abuik/i* — lO^"*). 



(ii) Exact diagonalizations (see section IV) that were limited to small sizes and large values of Abuik/i* (Abuik/-D = 
0.2). 
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FIG. 13. Self consistent gap profile with L = 128 sites and two values of io: to — 0.0625 (solid line) and to = 0.125 (dotted 
line). The other parameters are the same as on Fig. |l2|. The difference between the parallel and the antiparallel superconducting 
order parameters cannot be distinguished on the scale of the figure. 

In all approaches we find that the F/S/F trilaycr with metallic ferromagnets is characterized by Ap > Aaf while the 
F/S/F trilayer with insulating ferromagnets is characterized by Aaf > Ap. 

Finally, we mention two recent theoretical articles |^9[ in which Usadel equations have been used to treat F/S/F 
heterostructures in the diffusive regime. These authors have found that the metallic F/S/F heterostructure was 
controlled by pair breaking (Ap < Aaf) while we have found here an opposite behavior (Ap > Aaf)- In fact we 
believe that both approaches are correct but do not incorporate the same ingredients. The behavior of the model that 
we consider here is strongly related to non local pair correlations and can be explained with lowest order diagrams. 
The existence of a simple explanation shows the validity of the picture proposed in our article. On the other hand the 
model considered in Ref. p9[ is well suited to describe diffusive heterostructures. We are not certain at the present 
stage that the usual form of Usadel equations can be used to describe non local processes as we did here in ballistic 
systems. We think that a lot of understanding can be gained by discussing non local Usadel equations. 
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APPENDIX A: EXPRESSION OF THE KELDYSH PROPAGATORS 



The goal of this appendix is to rederive the main results of this article with the non equilibrium form ( |ll| ) of the 
Gorkov function, rather than the equilibrium Gorkov function (Ra). This formalism based on non equilibrium Green's 
functions is more general than the equilibrium Green's function formalism because it can also be applied to non 
equilibrium problems. A detailed investigation of this issue will be presented in the future. Here, we want to show 
that both formalisms coincide for the equilibrium problem, which constitutes also a test of the calculations presented 
in the main body of the article. 



1. One-channel problem 

Let us first consider the single channel model (see Fig. 0). The Green's functions G^'^ are the solution of the 
Dyson equation 



G° 
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FIG. 14. Variation of 5^ defined by (p4) with L = 128 sites in the superconductor. The parameters are identical to Fig. Iisl 



from what we deduce G"'" = gl'i/T), with T) given by (|I9). The Dyson-Keldysh equation associated to an arbitrary 
site in the superconductor takes the form 

^/5,,3 ^ 3pjj + 9p'.a ^a.aG^p + G f}^a^a,aga',p + G^gta^agala ta,aG^fj + Gp^^ta^aCJa^a ^a,aG^fj. (Al) 

Evaluating the five terms in ( [Al| ) leads to 
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Pf 

^(3.a,R fa,f}.A a,a 

y J Pi,i- 



The final step is to show that with fia = fJ-s this expression coincides with (|l8|). 



2. T^vo-channel problem with antiparallel magnetizations 

The Dyson-Keldysh equation associated to an arbitrary site /3 in the superconductor is the following: 

^/3,'/3 — 9p]l3 + 9p]a ^a.aGg^p + 5/3,'a'*a',a' ^a',/3 + G fj^J'a.aga.p + G'^,a'*a',Q'3Q,''^^ + G f^g^ta^ada^a ta,aG^f} 
+ Gi3a'ia',a'gJ^a'^a',a'Gg^,p + G /jg^ta^aQ a^a'^a' ,a'G ^, jj + G p^^ita' ,a' 9 a' ,of'^^°-'~^ ^^^0 
+ Gp^^ta^aiia.'a ^a.aG^jj + G f^^^.ta' y Cj ^,^^ita' ,a'G ^, jj. 

We need the expression of the following Green's functions: 



'^ — ' 9i,i 



QP.a' ^ _^a',a'^a> 



^/3,a Q 
jP.a Q 

/'''"' 

gP'"^' 



G"'^ = 



^2,1 "-^2,2 



(A3) 



(A4) 
(A5) 
(A6) 
(A7) 
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^a ,j3 



'-'l.l '-^1.2 



/ 



q',/3 ~a',/3 
•-^2,1 i/ 



(A8) 
(A9) 



We deduce from (e3) - (EQ) and ( A4) - (A9) the final form of the Gorkov function in the antiparallel ahgnment: 
^/3,/3 - ^«7i^'^-f('^ - Msj |P/ + K I 5ia / Pg + K I 52,2 3 P/ + K I 5i,i 5 P/ (AlUj 

I ua',Q' |2 a'.a',-R fl3,a',Ra'.p , l .a,Q 1 4 a,a,i?a,a, A --/3,q._R ?a./3, A Q.a , l .a'.a' |4 a'.a',-R a',o',A~Q',/3.A fP.a'.Ra'^a' 

"T I*- I i/2,2 i Pg "T r I i/i,i yi,i y i "g ^ \ I ^2,2 1/2,2 9 J Pg 

I I j.a,a |2 I j.a',Q' 1 2 a,a,fla',a', A p-/3,Q,fl,-.Q',;3, A a,a' , l .a,a |2 l .a', a' |2 a,a,Aa',a',i? ?/3,Q',fl ?q,/3,A «',« \ 

"'"I'' II'' I yi,i y2,2 y y P/ "•" K I r I yi,i 1/2,2 j j ^/ j 



Using the propagators obtained in section [II B 1 to evaluate ( A10| ), we can show that Eq. ( A1C| ) leads directly to ( p^ ) 
and d^c" 



3. Two-channel problem with parallel magnetizations 

Let us now consider two single-channel ferromagnetic electrodes having a parallel spin orientation. Following 



section 



A 2 , we obtain 



^/3,/3 - 2?7r7T,_F(w-/isj |p^ + I* I ffl,l / Pg +K I 5l,l / Pg + I* I 5l,l 5 P/ (Allj 

I ua',a'|2 a',a',_R~;3,a',fl a',/3 , l .a,a |4 a,a,fl a,a,Ap-/3,a,i?, jq,/3,A Q,a , l j.a',Q' |4 a',a',i? a',a',A ?a',;3,A ~/3,q',_R «',»' 

+ I' I yi,i y P/ "1" I' I yi,i yi,i y / Pg ^ \'' I yi,i yi,i J y Pg 

1^ |j-a,a|2|j.a',Q' |2 a,a,-R a',a',A-./3,a,i?, ?q',/3,A a'.a i^ l ia,a |2 1 .a'.a' |2 a,a,A a',a',_R--/3,Q',fl ?q,/3,A q',q\ 

"T I*- I r I yi,i yi,i y J Pg "•" I' I r I yi,i yi,i y J ^a \ 

+ 2tTrnF(i^ - Pa)\t \ g J Pij^ + 2nrnp(uj - ^a')\t \ 9 J Pi,i , 

where the propagators g and / are given by 

j-a,0 ^ |^a',a'|2^aV /a,a'j-a'.,0 _ ^a' .a' j-a.A 



r"- 



Vy 



fa,0 ^ J_ 



/■" = 



Vf 



■' Vf 



B.a I ua', q' |2 a .a ( a'.a 3, a' a'.a' 0,a\ 

9 +\t ' I 52,2 [9 9 -9 9 ) 



where I?f is given by Eq. (B2). We can show that Eq. (All) leads directly to (pO|) and (& 



(A12) 
(A13) 
(AM) 
(A15) 



APPENDIX B: FERROMAGNETIC ELECTRODES WITH AN ARBITRARY SPIN POLARIZATION 



1. Dyson matrix 

The advanced and retarded Green's function Gp^p at an arbitrary site f3 of the superconductor can be deduced from 



Eq. 



/^/3,/3 _ „/3,/3 , ^P,a,a,ar~ia,(3 B.a.a.ar^a.p . 0,a' ,a' ,a' ^a' ,13 d.a' ,a' .a' ^a' ,f3 

"-^2,1 —52,1 +52,1'' '-'1,1 ~ 92,2^ "-^2,1 +52,1 '' '-'1,1 ~ 92,2 ^ '-'2,1 ■ 



(Bl) 



Eq. (Bl) can be used to evaluate the equilibrium Gorkov function given by Eq. (O) and deduce the value of the 

self-consistent superconducting order parameter. The Green's functions G^ i , Gj']^, G" '{ and Gj { are the solution 
of the Dyson equation (pi) which can be expressed as a 4 x 4 Dyson matrix: 
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1 - K^'"t"'' 

1\2,1 t 
■"^2.1 '' 



K^ft° 



'"1,2'' 



1 - K'^-ft^'" 



Kl^tP-'' 




r r<°;ii 1 
<^i,i 




-"■1,1 


-K^itP^^ 




<-^2,l 




--f*^2.1 


K\itP^^ 




<^i,i 




-"■1,1 


- Klitp-^ 




<^2,l 




-"■2,1 



(B2) 



The coefScients Kj j are of the form iiTij = gi^itgij. For instance, -?ir°'" = glf ^""'^ gtT ■ The inversion of Eg. (|B2| ) is 
described in section B2 for f-wo ferromagnetic electrodes having a parallel spin orientation, and in section B2 for two 
ferromagnetic electrodes having an antiparallel spin orientation. 



2. Parallel magnetizations 



If t'*'" = t^'P ^ the 4x4 Dyson matrix given by Eq. (B2) can be -written in terms of 2 x 2 blocs: 

L*p Kp 



(B3) 



■where Kp and Lp can be obtained from Eq. (B2): 

l-i^l,l i^l,2 

^ K2.1 1 - K2a 



and Xaf 



— ii,i -^1,2 
^2,1 —-^2,2 



The inverse of 2?f given by (B3) takes the form 



t)-^ - 
Up — 



-Kp^L*pMp^ 



-Kp^LpNp' 



■with Mp = Kp — LpK-p ^Lp and Np = Kp — LpK^ ^Lp. The matrix in Eq. (B3) can be evaluated explicitly to obtain 
the different terms in the Green's function (Bl): 



1 



((Im g^:rt"''^Gl:', ))^-^nplf{A^,0) 



x^ 



ao 



((Im -52^"i"'''G^f )) = - 



■a,/3 



(1-A'i,i)2(l-X2,2)2 



1+2^1 \Raj3J {^ + X]Y Ra^pRa' .pRa,, 



9^:2 {-^2,1(1 - i^l,l)(l - i^2,2)i^tf 



(B4) 
(B5) 



(1 - K2,2Kl - Ki,ifk^'f - [k2,i (ii,i(l - ^2,2) + ^2,2(1 - i^l,l) 



1^ 



i2,l(l-i^l,l)(l-i^2,2) 



K 



b.l3 



l~ 



^2 r>b.P 



^^^--L2.2{l~K,.,YKl'^, 



((Im [gP'^ft^^'^^'Gif 
((Im \~g'2at'''''''Gif 



» 



:^Po/(Ac.',/3) 



X] 



l+X^ \Ra'Ji 



flo 



(1 + x'i Ra,fiRa',f5Ra,a' 



{l-Kisf{\-K2,2? 



f. 



"' { [-^2,1 (il,l(l - -?^2,2) + -£2,2(1 - Xi,i) 



(B6) 

(B7) 



1 ^ 



^2,1(1 --?^l,l)(l--?^2 



K 



a, 13 



1 ~ 



\2i>a.,p 



la - 2^2'2(1 - K.^iYK^^ 



~ i^2,i(l - ifi,i)(l - ^^2,2)-^-!,? - ^(1 - -?f2,2)(l - K.^ifklf 

■where X| and x^ are given by Eqs. (R) and (ph. The function /(A) is given by /(A) = /S./^/lo'^ — A^. The phase 
contribution has been factored out in the coefficients Lij. For instance, Li,i — Li_i exp (i(/3Q,^) exp (i^'a./s), 'with 



Li^i = -Ti'^t^p^ p^ {aa/Ra.p)- To obtain Eq. (^) -we use a simplification of Eqs. (^4|) - ( p37| ) in -which Eqs. ( p34| ) 
(B7) are transformed into a local equation. This is done by replacing A^^ and A^/^ ■with A^. 
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3. Antiparallel magnetizations 



The case of an antiparallel spin orientation of the ferromagnetic electrodes can be treated in a similar manner. 
Once the lines and columns 3 and 4 have been interchanged, the Dyson matrix given by Eq. (B2) takes the form 



with 



The inverse is given by 



Kaf = 



2?AF 



1 - Ki,i Ki,2 

K2.I 1 - if 2,2 



-'^AF ^AF 



and Laf = 



Li,2 ^Lii 

— i2,2 -^2,1 



-^AF 



Maf 



^k^^lmmXp)-' 



-(^AF)-'iAF^AF (MIj.)-' 



with Maf = ^AF — LApiK'p^-p) ^L\p. The different terms in the Green's function (Bl) are the following: 

"2 



1 



((Im g!^:rt"'^Gl:^ )) = --np'omo.,,) 



x^ 



ao 



^T^i 



"0 



((Im -52T*"'"G^f » = - 



l + Xj \Ra.j3J {1 + X^){1 + Xl) Ra,pRa',f3Ra,a' 



(l-Xi,i)(l-if2,2) ' 1 ' 2^ ' ' 2'^ 2 ' ^'1 



1 



„P-a' ,a' ,a' /^a ,[j 



((Im g,yr'--Gti'^ )) = -^vrpo^/(A„,,^) 



1 - if2,2 
1 



^^2, 1-^1,1 + Ki^2L2.2 + ^^24(1 — -fi'l.ly 



if 



b,/3 



Oo 



l+Xi \Ra',[i) {l+X'\){l+Xi) Ra.f3Ra',l3Ra.a' 



0,a' ,a' .a' /^a' ,13 



Ki,2L2,2 + K2,lLi^i + -Li,2{l — K2,2) 



+ iii,i(l - ifi,i)(l - K2,2)k^2f + ^(1 - ^l.l)(l - ^2,2)'if2',f 
+ ifi,2(l-ifl,l)(l-if2,2)^l.f}. 

Using a "local" approximation in which A^,^ and A^/^ are replaced by A^ leads to Eq. (15. 



(B8) 
(B9) 

(BIO) 
Kll^ (Bll) 



-a,/3 
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